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Abstract 

Quantum branching programs (quantum binary 
decision diagrams, respectively) are a convenient 
tool for examining quantum computations using 
only a logarithmic amount of space. Recently sev- 
eral types of restricted quantum branching pro- 
grams have been considered, e. g. read-once 
quantum branching programs. This paper con- 
siders quantum ordered binary decision diagrams 
(QOBDDs) and answers the question: How docs 
the computational power of QOBDDs increase, if 
we allow repeated tests. Additionally it is de- 
scribed how to synthesize QOBDDs according to 
Boolean operations. 

Keywords: Computational Complexity, The- 
ory of computation. Quantum computing. Branch- 
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1 Introduction 

A central question of quantum complexity is, in 
which cases quantum computations do outper- 
form classical ones. Famous results are the al- 
gorithms of Shor (Cni) and Grover ([7]); apart 
from that much has been achieved by examining 
various restricted models of quantum computation 
and comparing them with their classical counter- 
parts. One such model are branching programs 
(BPs). They are related to circuits. Boolean for- 
mulas, and nonuniform space complexity. Quan- 
tum branching programs have been considered in 

isi, la, in, 0, in and 0. 

A deterministic BP B on the variable set 
{xi,X2, ■ ■ ■ jXn} consists of a directed acyclic 
graph G = (V, E) where multi-edges are allowed. 
Two of the nodes are denoted as targets. They are 



sinks in the graph-theoretical sense, and are la- 
beled and 1. The other nodes are called branch- 
ing nodes. They get labels from {xi, X2, . . . , Xn}. 
The edges get labels from {0, 1}. For each branch- 
ing node, there is exactly one outgoing edge la- 
beled 0, and one outgoing edge labeled 1. The 
size of a BP is the number of its edges. 

A branching program computes a Boolean func- 
tion in a natural way. Each input a G {0, 1}" 
activates all a^-edges leaving aii-nodes, for i = 
1, 2, ... 71. A path in G is defined to be activated 
by a, if a activates all its edges. An input a is 
accepted if the path activated by a leads to the 1- 
sink and is rejected in the other case. This model 
can be generalized to nondeterministic or proba- 
bilistic modes of computation in a straightforward 
way, see [T^ . 

Quantum branching programs (QBPs) can be 
defined by adding transition amplitudes to the 
edges and allowing more than two sinks (see [HI). 
We outline this approach very briefly. The com- 
putation on an input a starts at the source of the 
QBP. With respect to the transition amplitudes 
each step of the computation consists of a super- 
position of nodes. Finally a measurement deter- 
mines the result of the computation, i.e. the la- 
bel of the resulting sink. Certainly the transition 
amplitudes have to fulfill a global well-formedness 
constraint that ensures a unitary evolution of the 
computation. 

In this paper we consider another - equivalent - 
approach following 0. This approach is particu- 
larly useful for leveled branching programs, where 
the nodes are partioned into levels Li, . . . , Lm. Lm 
consists of the sinks, for i < m, the nodes in Li 
are labeled by the same variable and the outgoing 
edges of a node in Li lead to nodes in the level 
below, i.e. in L^+i. 

A quantum branching program on the variables 



{xi, X2, ■ ■ ■ , Xn} of width w and length I consists 
of 

• a set D of cardinahty w; we assume D = 

m,...,\w)}, 

• a sequence of pairs of unitary tranforniations 
{Ty. , Ty. ),i = 1, . . . ,1 on the complex vector 
space spanned by D where j/i is a variable in 

{xi,X2,-.-,Xn} 

• a starting state \1) £ D and 

• a set of accepting states F C D. 

We call the sequence {yi, . . . ,yi) of variables in 
{xi,X2, . ■ . ,Xn} the variable ordering of the QBP. 
The computation proceeds in the complex vector 
space spanned by D. Each state of the computa- 
tion is a vector of length 1 in this space. The 
computation starts with state On an input 
a £ {0, 1}" the transformation T^^ is applied if 
a assigns the variable yi to 0, in the other case 
Ty^ gets used. The result is a unit vector on that 
either Ty^ or Ty^ is applied. After I steps the com- 
putation stops by measuring the state 

where a assigns yi with Ci, and {ipi^a)) is a vec- 
tor (ai, whose components are complex 
amplitudes, or - equivalently - the superposition 
127=1 Q^il*)- The measurement results with proba- 
bility |aip in the state \i). If this result is a mem- 
ber of F we accept the input, in the other case 
we reject. Therefore the state \ipi{a)) plays an im- 
portant role - this leads to the definition of the 
final amplitude in a computation. Let B be the 
QOBDD defined above. Then the final amplitude 
of \i) according to the input a is 

finalAmp(|i), a) :— {i\ilJi{a)), 

i.e. the component of |i) just before the measure- 
ment; (i I j) denotes the inner product of the com- 
plex vectors \i)\j). For each \i) £ D the measure- 
ment finishing the computation of S on a yields 
the result \i) with probability | finalAmp(|i), a)p. 
An input a is accepted with probability 

^|finalAmp(|i),a)p. 

Sauerhofi^ and Sieling proved in that QBPs of 
polynomial size correspond to logarithmic space 
restricted computations of nonuniform quantum 
Turing machines. In our model the size is the 



product of width and length. Ablayev, Moore and 
PoUett proved that NC^ can be accepted by QBPs 
of width 2 and polynomial length, see Upper 
and lower bounds have been proved for several re- 
stricted versions. An important variant are quan- 
tum ordered binary decision diagrams (QOBDDs). 

Branching programs are important not only in 
theory but also in applications. In this con- 
text they are denoted as binary decision diagrams 
(BDDs). BDD-based data structures for Boolean 
functions play a key role in hardware verification, 
test pattern generation, symbolic simulation, logi- 
cal synthesis or analysis, and design of circuits and 
automata (for a survey see 22). Once the model is 
chosen one needs efficient algorithms for many op- 
erations, particularly for synthesis, minimization 
and equivalence test. The non-equivalence test for 
two functions / and g is equivalent to the satisfia- 
bility problem for f ® g. It is known that the sat- 
isfiability problem for read-twice branching pro- 
grams is NP-complete. Therefore one prefers the 
restricted types of branching programs. Very im- 
portant is one introduced by Bryant |H] that may 
be regarded as the state-of-the-art data structure 
in many applications. 

A QBP as defined above is a QOBDD if the 
variable ordering (j/i, . . . ,yi) is a permutation of 
{xi, X2, ■ ■ ■ , Xn}- Or, more illustrative, the length 
is n and each uniform transformation depends on 
another variable. 

QOBDDs have been considered by Sauerhoff 
and Sieling in |5] . They have presented a function 
computable by succinct QOBDDs that requires 
exponential size deterministic OBDDs. Counter- 
wise they have found a very simple function that 
is not computable by polynomial size QOBDDs. 
They call this function N0„ (neighbored ones). 
It is defined on n variables {xi, a;2, . . . , and 
tests whether there are neighbored variables with 
value 1, i.e. an input is accepted if and only if 
Xi = Xi^i ~ 1 for some i < n. This function is 
computable by deterministic OBDDs of size O {n). 
This weakness of QOBDDs has the reason that ev- 
ery step of the computation is a unitary and there- 
fore reversible transition. For strongly restricted 
models of quantum computations (the situation is 
similar for some kinds of quantum finite automata) 
it seems to be difficult to forget variables already 
read (see PI). 

Thus the question arises: How does this situa- 
tion change if we slightly diminish the restriction? 
How does the computational power change if we 
consider QOBDDs with repeated tests? This leads 



to the definition of k-QOBDDs. 

Unformally a fc-QOBDD is the concatenation 
of k QOBDDs that obey the same variable or- 
dering. More presicely a fc-QOBDD B accord- 
ing to the variable ordering cr is a QBP B on 
the variables {xi, . . . , Xn} of length kn. The lev- 
els are partitioned into k layers of length n. In 
each layer the variables are tested according to 
a. The computation of a fc-OBDD on an input a 
is determined by a sequence of unitary transfor- 
mations Tl, . . . , T^„, e e {0, 1}, where T^^ is cho- 
sen according to the input bit a(j(i)- Thus, on 
input a the i-th layer performs the computation 

rpO-^(n) rpO-<T(2) rp<i^(l) 

BoUig, Sauerhoff, Sieling and Wegener have 
proved in ^ a hierarchy on deterministic fc- 
OBDDs. It turns out that the computational 
power of polynomial size fc-OBDDs is strictly 
greater than that of (fc — l)-OBDDs. Things are 
different in the case of nondctcrministic fc-OBDDs. 
In |5] it is shown that nondctcrministic , parity 
and randomized fc-OBDDs are not more powerful 
than OBDDs with the correspondent computation 
modes. 

Section [21 starts with an analysis of the way a 
fc-QOBDD computation on an input a evolves. To 
apply this analysis for comparing QOBDDs with 
fc-QOBDDs, we consider products of QOBDDs in 
13] This provides a method to perform the Boolean 
synthesis. Scction^shows that repeated tests are 
of no use in the case of QOBDDs with unbounded 
error. 

2 The way a A;-QOBDD com- 
putes 

We consider a fc-QOBDD B with variable order- 
ing a and width it; on n variables. The computa- 
tion evolves according to the unitary transforma- 
tions T{, . . . , Tl^, e e {0, 1}. We define U^{a) to 
be the transformation performed by the i-th layer 
under a. Formally, U,{a) = T^;'"' • . . ■ • T(":'j;.„+2 • 
T^i'^iyn+i- The final amplitudes of the computa- 
tion on a are the components of the superposition 
Ukia)- ...■Ui{a)\l). 

Let a[j^ (a) be the amplitude of \j) in the state 
Ux{a)\i), i.e. 



functions from {0, 1}" to C by 

t,f\a):=a['J{a), j e {1,2, . . . ,w}. (1) 

For As {2, 3, ... , fc}, let fi'^^'' be a u; x w-matrix 
of functions from {0, 1}" to C defined by 

^J!^(a}:^Q^^{a), j & {1,2, . . . ,w}. (2) 

According to our definition, ^['^j(a) equals the 
amplitude of \j) in the result of the computation 
of the A-th layer on a starting with 

For |i) E D we define f3i := finalAmp(a, An 
easy calculation reveals that for all a G {0,1}" the 
vector of final amplitudes Pi can representated as 
a matrix product: 

(/3i(a),...,/3,„(a))^ = 
fi^^\af ■ fi^^\a) ■ . . . ■ ^l^''-^\a) ■ n^>'\a) (3) 

Figuring out the right hand side of Equation 13 
we obtain for j ~ 1, . . . ,w 

PM) = 

Z2,...,ifce{i,...,ii>}"' ~~ ' 

(4) 

We define the acceptance probability of B on 
some input a by 

acc{a) := ^ l/^Jp- 

Our purpose is to construct a quantum OBDD B' 
that simulates the quantum fc-OBDD B in the case 
of bounded error computations. B^ accepts an in- 
put a with probability greater than 1/2 if and only 
if acc{a) > 1/2. To this end we adopt the well- 
known "product-construction" for finite automata 
common for synthesizing BPs. 

3 Product construction and 
Boolean synthesis 

In the quantum case it is convenient to use the ten- 
sor product. Let Bi and B2 be quantum OBDDs 
using the transformations (T!j?. , T^.), i = 1, . . . ,n 
and {S^., S^.), i = 1, . . . ,n, respectively (the same 
approach works for fc-QOBDDs). Bi is defined on 
the set Di of cardinality Wi, i = 1,2. We denote 



We define the column vector /i*^^^ of length w of 



Bi®B2 



to compute on the set Di x D2 of elements \i) (E) \ j) 
by the transformations 

(r° ®5°,,T,i_®^ij, t = l,...,n. 

Let \di) be contained in D.i,i = 1,2. The com- 
plex values finalAmp(|c?i), a) are the according fi- 
nal amplitudes of the computations oi Bi,i = 1,2. 

Then for the computation of Bi ® B2 it holds 
that 

finalAmp(|(ii) ® 1^2), a) = 

finalAmp(|o?i), a) • final Amp ( 1^2), a). (5) 

Using Equation [S] and standard techniques 
as OBDD-probability-amplification we can prove 
that the logic synthesis operation is feasible for 
QOBDDs with different error bounds. 

Proposition 1 Let fi,i = 1,2 he functions com- 
putable by quantum OBDDs of width Wi. Then 
/1A/2 is computable by a quantum OBDD of width 
polynomial in wiW2- 

The proof of Proposition ^ is straightforward. 
We define the accepting states as Fi (X) F2. 
If Bi accepts an input a with probability pi,i = 
1, 2 then Bi ® B2 accepts a with probability piP2- 
Thus synthesizing two QOBDDs with error bound 
e result in a QOBDD computing the conjunction of 
the input QOBDDs with error bound 1 — (1 — e)^. 
By a finite number of additional synthesis steps 
the error can be decreased to e. 

4 Quantum fc-OBDDs with 
unbounded error 

We make use of the product construction de- 
scribed in the preceding subsection to simulate 
a fc-QOBDD by a QOBDD. Consider the fc-wise 
product of = . . . , \w)}, i.e. 

fc 

^® (g) ^ = {|h*2 ■■■Ik); \ii) eD,l^i,...,k}. 

i=l 

\i1i2 ■ ■ - ik) is the common abbreviation of (g) 
. . . ® \ik). To simulate the fc-QOBDD B wc de- 
fine a QOBDD B"^ computing on the set D"^. Its 
transformations are (Tf Tf '^), . . . , (T®-", T®'^) 
where T^''^,e G {0,1} are chosen according to 
Xa{n) ■ We define 

Tf''=Ti^T,^_,^^...^T:^^,_,),^, (6) 



for e G {0, 1}, i = 1, . . . , n. Note that on some in- 
put a the QOBDD performs the unitary trans- 
formation C/®(a) = C/i (a) (g) . . . (g) Uk{a). 

Let us examine how B"^ simulates the way B 
computes. Let a S {0, 1}" be fixed. We apply 
U^{a) on \li2i3 ■ ■ ■ ik) , where 12*3 ■■■ifc arc arbi- 
trarily chosen elements of D. Let 

'4'li2i3---ik i'^} • 

C/®(a)|1^2^3...^fc) = 
Ui{a)\l) ® U2{a)\i2) Uk{a)\ik) (7) 

We start the computation in state 1 11213 ■ ■ - ik) 
as above, then the component |«2*3---ifei) for 
\j) e -D of the state il^n^i^^Akio) has the same 
amplitude as thefoUowing computation path tt of 
B. TT starts with |1), the intermediate result after 
the first layer is ^2, after the second layer we reach 
^3 etc; the final result of the considered path tt is 
Note, that it is quite natural to carry the con- 
cept of a superposition of states to a superposition 
of computation paths. Formally it holds that 

(1213 . ..iki I V'H2»3...JA,(a)) = Aii2«...ifci(a): (8) 

using the notation of Equation 0] 

We utilise Equation|Hlto build a QOBDD 6' that 
performs the same computation as the fc-QOBDD 
B. We define D' := U {\to),\ti)}. ti is an 
accepting state and to rejecting (in the following 
we often abbreviate \j) G Z? by j G D.) 

F' := {\i2i3---ikj);i2,i3,---,ik S D,j G F} 

is the set of accepting states of B\ Let V he a. 
unitary transformation from D' on itself that maps 
the vector |1 . . . 1) to the superposition (let m := 

V\l...l) = 

= ^ E |l*2,...,ifc) +:r^|to) 

V2m . 2y/rn 



The images of all other members of Z?® are cho- 
sen such that V is unitary. This is possible, since 
y|l . . . 1) is a vector of length 1. 

Let i G {1, . . . , n}, e G {0, 1}. The transforma- 
tions T®'*^ on are defined according to Equa- 
tion |5| We define the unitary transformations T^'^ 
from D' to itself as behaving on as J"®'^ and on 



{to,ti} as the identity: T-^ld) =Tf\d) for d e D® , 
T;^|to) = |io)and7;'^|ti) = |ti). 

Now is defined as computing according to 

{V -t'.W ■t[%{t'^\t'^%...,{tI^X')- 

We determine the acceptance probabiUty of the 
computation of B' on an input a. i3' starts with 
state |1 . . . 1). Applying V on this start state has 
the result described in Equation]^ Thus, applying 
U'{a) := T^<"' 
has the result 



.T^'^^Kt^'W -Fon 



1 



'2m 



J2 U'{a)\U2,..., 



»2 ik€D 



1 , , \/2m — 1 , ^ 
+ :^^|^o)+ „ ^ \ti). 



The first part of this sum can be rewritten as 
1 

i2,...,ikeD jeD 

1 



/2m 




J2 



12 ifce-D 



+ 



For the acceptance probability this implies 
accg' (a) = 

-E 

12 jfceJ 



V2m- 1 



2V?^ 
2to- 1 



2m 



4m 



where (3j{d) is defined as in Equation 0] Observe 
that the acceptance probability of the fc-QOBDD 
S on a is Y.j(^F 

Thus, if B accepts a with probability at least 
1/2 then S' will do the same. We formulate this 
result as 

Theorem 1 For all natural numbers k > I a 
sequence of Boolean functions {fn) computable 
by polynomial size quantum k-OBDDs with un- 
bounded error is also computable by polynomial 
size quantum OBDDs. 

References 

[1] F. Ablayev. On comparative computational 
power of classical and quantum branching 
programs. Technical report, BMPI, 2003. 



[2] F. Ablayev, A. Gainutdinova, and M. Karpin- 
ski. On computational power of quantum 
branching programs. Proceedings FCT 2001, 
LNCS 2138:59-70, 2001. 

[3] F. Ablayev, C. Moore, and C. PoUett. 
Quantum and stochastic branching pro- 
grams of bounded width. Technical report, 
arXiv:quant-ph/0201139, 2002. 

[4] B. BoUig, M. Saucrhoff, D. Siehng, and 
I. Wegener. Hierarchy theorems of fcOBDDs 
and /cIBDDs. Theoretical Computer Science, 
205:45-60, 1998. 

[5] H. Brosenne, M. Homeister, and St. Waack. 
Nondcterministic ordered binary decision dia- 
grams with repeated tests and various modes 
of acceptance. Information Processing Let- 
ters, in press, 2004. 

[6] R. E. Bryant. Graph-based algorithms for 
Boolean function manipulation. IEEE Trans- 
actions on Computers, 35:677-691, 1986. 

[7] L. Grovcr. A fast quantum-mechanical algo- 
rithm for database search. 28th STOC, pages 
212-219, 1996. , arXiv:quant-ph/9605043. 

[8] M. Nakanishi, K Hamaguchi, and T. Kashi- 
wabara. Ordered quantum branching pro- 
grams and more powerful than ordered proba- 
bilistic branching programs under a bounded 
width restriction. 6th COCOON, LNCS 
1858:467-476, 2000. 

[9] M. Saucrhoff and D. Sieling. Quan- 
tum branching programs and space-bounded 
nonuniform quantum complexity. Theoretical 
Computer Science (to appear), arXiv:quant- 
ph/0403164, 2004. 

[10] P. Shor. Polynomial-time algorithms for 
prime factorization and discrete logarithms. 
SIAM Journal on Computing, 26:1484-1509, 
1997. arXiv:quant-ph/9508027. 

[11] R. Spalek. Space complexity of quantum com- 
putation. Master's thesis, Karl's University 
Prague, 2002. 

[12] I. Wegener. Branching Programs and Binary 
Decision Diagrams - Theory and Applica- 
tions. SIAM Monographs on Discrete Math- 
ematics and Applications. SIAM, Philadel- 
phia, 2000. 



